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Abstract 

We introduce the concepts of complex Grassmannian codes and designs. Let 
Gm,n denote the set of m-dimensional subspaces of C": then a code is a finite 
subset of Qm,n in which few distances occur, while a design is a finite subset 
of 

Qm,n that polynomially approximates the entire set. Using Delsarte's linear 
programming techniques, we find upper bounds for the size of a code and lower 
bounds for the size of a design, and we show that association schemes can occur 
when the bounds are tight. These results are motivated by the bounds for real 
subspaces recently found by Bachoc, Coulangeon and Nebe, and the bounds 
generalize those of Delsarte, Goethals and Seidel for codes and designs on the 
complex unit sphere. 

1 Introduction 

In this paper, we introduce the concept of complex Grassmannian codes and 
designs: codes and designs in the collection of fixed-rank subspaces of a complex 
vector space. 

In the 1970's, Delsarte [10] developed a series of excellent bounds for certain 
error-correcting codes by treating codewords as points in an association scheme 
and then applying linear programming. Shortly thereafter, Delsarte, Goethals 
and Seidel [11] showed that the same technique could also be used on systems 
of points on the real or complex unit sphere, which they called spherical codes 
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and spherical designs; this resulted in important contributions to problems in 
sphere-packing [9l Chapter 9]. This linear programming technique, which is 
now known as "Delsarte LP theory", has proved surprisingly portable. Re- 
cently, Bachoc, Coulangeon and Nebe [3] generalized the results of Delsarte, 
Goethals and Seidel to real Grassmannian spaces, and Bachoc [2j pointed out 
that "the same game" can be played over the complex numbers. In this paper, 
we investigate more closely the case of complex Grassmannian codes. 

The motivation for studying complex Grassmannians comes from the the- 
ory of quantum measurements. Roughly speaking, any complex Grassmannian 
1-design defines a projective measurement in the theory of quantum mechanics. 
It has recently been discovered that complex projective 2-designs correspond 
to quantum measurements that are optimal for the purposes of nonadaptive 
quantum state tomography [21j. In fact, this is also true in the more gen- 
eral Grassmannian setting: complex Grassmannian 2-designs are the optimal 
choices of measurements for nonadaptive quantum state tomography when the 
observer only has access to measurements with a restricted number of outcomes. 
More details will appear in a paper by Godsil, Rotteler, and the author |13j . 
Complex Grassmannians also play a role in certain wireless communication 
protocols [1]. 

Define Qm,n to be the set of m-dimensional subspaces of an n-dimensional 
complex vector space. Without loss of generality, we will always assume m < 
n/2. Usually, we will represent a subspace a by its n x n projection matrix Pa- 
The inner product on Qm,n is th.6 tracG inner product for projection matrices: 

{a,b) :=tr(P;n) 

Since (a, 6) = {b,a), the inner product is real. This is a measure of separation, 
or distance, between two subspaces — note that is not a distance metric per se: 
the inner product of Pa with itself is maximal rather than minimal. However, 
the chordal distance [8], defined by 

dciPa,Pb) := Vm-ti{PaPb), 

is a monotonic function of the inner product. Given a finite set of inner product 
values A, an A-code is a subset S of Gm,n such that 

A = {tv{PaPb) ■■a,b£ S, a^b}. 

An s-distance set is an ^-code with |^| = s. This generalizes the concept of 
an s-distance set on the complex unit sphere: if u and v are unit vectors, then 
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their separation distance on the unit sphere is a function of 

I * 1 2 , / * * \ 

\u v\ = tr(nn vv ). 

We are interested in codes of maximal size for a fixed ^ or s, and bounds on 
their size based on zonal polynomials. Table [T] in Section [U] gives a summary of 
the bounds for small |^|. 

The outline of this paper is as follows. In Section [21 we describe the orbits 
of pairs of subspaces in Qm.n under the action of JJ{ti)'. these orbits play a 
significant role in the bounds derived later on. In SectionsEHHandEl we develop 
the necessary representation theory background needed for our LP bounds. 
In particular, we discuss the decomposition of the square- integrable functions 
on Qm,n into irreducible representations of U{d), and the zonal polynomials 
for these representations. The results in this section are all known, and the 
development is quite similar to that of Bachoc, Coulangeon and Nebe for real 
Grassmannians. In fact, the complex case is actually easier than the real case, 
because representations of the unitary group U (n) are easier to describe than 
representations of the orthogonal group 0{n). In Section[6l we develop absolute 
and relative bounds for codes, and show how these bounds for Qm,n reduce to 
known bounds for complex spherical codes when m = 1. These bounds are 
compared to some other known bounds for subspaces in Section [71 In Section [H 
we consider Grassmanian designs. Grassmannian codes enjoy a form of duality 
with complex Grassmannian designs, very similar to real Grassmannian codes 
or spherical codes. In Section [HI we give examples in which the bounds are 
tight. In many cases codes of maximal size or designs of minimal size have the 
structure of an association scheme, which we describe in Section [TOl 

2 Orbitals 

In this section we describe the orbits of pairs of elements of Qm,n under the 
action of U{n). 

First, we claim that Qra,n can be identified with a factor group of the unitary 
group, U{n) /{U{m) x U{n — m)). For, consider the first m columns of a matrix 
of U{n) as the basis for a subspace a of dimension m in C", letting the last 
n — m columns be a basis for a^. Then a is invariant under the action of U{m) 
on the first m columns, while O"*" is invariant under U{n — m). 

As a result of this factor group, U{n) acts on Qra n follows: if U is in 
U{n) and Pa is the projection matrix for a £ Qm,n, then 

U-.Pa^UPaU*. 
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This action is an isometry, in that it preserves the trace inner product on Gm,n- 
Unhke the complex unit sphere, however, U{n) is not 2-homogeneous on Qm,n- 
U{n) does not act transitively on pairs of subspaces with the same distance. In 
other words, the fact that tv^PaPb) = ti{PcPd) does not imply that there is a 
unitary matrix mapping a to c and b to d. In order to use zonal polynomials, 
we need to understand the orbits of pairs in Gm,n under this isometry group, 
which requires principal angles. 

Given a and b in Gm,n, the principal angles 9i, ... ,9m between a and b are 
defined as follows: firstly, 9i is the largest angle that occurs between any two 
unit vectors ai £ a and 61 G b: 

61 := min arccos I at 61 1 . 

ai Ga 
bieb 

Secondly, 62 is the largest angle that occurs between any two unit vectors 
a2 G an a;/- and 62 G bnbj^. Similarly define 9^, ... ,9m- These principle angles 
are closely related to the eigenvalues of PaPb- the first m eigenvalues of PaPt are 
{cos^ 01, ... , cos^ 9m}- Because of this correspondence, for the remainder of this 
paper we simply refer to the eigenvalues yi := cos'^ 9i (rather than the values 
9i) as the principal angles between a and b. Note that n — moi the eigenvalues 
of PaPh are zero, so we need only consider the first m eigenvalues. Conway, 
Hardin, and Sloane [8] accredit the following lemma to Wong |241 Theorem 2]. 

Lemma 2.1. The principal angles characterize the orbits of pairs of subspaces 
under U (n). 

Proof. Suppose U € U{n) maps projection matrices Pa and Ph to Pc and Pd 
respectively. Then by similarity, the eigenvalues of 

PcPd = {UPaU*){UPbU*) = UPaPbU* 

are the same as the eigenvalues of PaPb- 

Conversely, we show that if PaPb and PcPd have the same eigenvalues, then 
some unitary matrix U maps a to c and b to d. We do this by unitarily mapping 
a and b into a canonical form that depends only on the eigenvalues of PaPb- 

Let Ma be an n X m matrix whose columns [ai, . . . , Um] are an orthonor- 
mal basis for a, so that MaM* = Pa and M*Ma = I- Similarly define 
^b = [bi,---,bm] for b. Suppose M*Mb has singular value decomposition 
UDV* , where U and V are m x m unitary and D is m x m diagonal. Then 
{MaUy* (MbV) = D. Since the columns of Mall are another orthonormal basis 
for a, without loss of generality we replace Ma by Mall and likewise replace 
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Mb with MfjV. In other words, we may assume without loss of generality that 
M*Mb = D, where D is a diagonal matrix of singular values. 

Next, define the columns of Na = [om+i, • • • , On] to be any orthonormal 
basis for a-*-, so that NaN* = I — Pa and N*Na = /. Further assume that 
N*Mj) = QR, where Q is {n — m) x {n — m) unitary and R is {n — m) x m 
upper triangular (the Qi?-decomposition of N*Mfj). Then Q*N*Mi, = R, and 
the columns of NaQ form another orthonormal basis for a^. Replacing Na by 
NaQ, we may assume without loss of generality that N*Mh is upper triangular. 

Finally, let Ua ■= ^ tv* ) ' ^^^^ is an n x n unitary matrix. Then 

UaMa=(^^^y, UaM,= (^^y 

If PaPb has eigenvalues cos^^j, then M*Mh = D has singular values cos^j. 
Moreover, since UaMi^ has orthonormal columns, it follows that R also has 
orthogonal columns. We may therefore assume that R is not just the upper 
triangular but diagonal, with diagonal entries sin^j. Thus Ua is a unitary 
matrix which maps Ma and Mj, into the form 

cos Oi \ 
cos 9m 

sin 9i 

sin 9m 
/ 

Since any pair {Ma,Mi,) with principal angles cos^ Oi can be mapped to this 
canonical form, it follows that the eigenvalues of PaPb characterize the orbits 
of pairs (a, 5) under the unitary group. □ 



3 Representations 

In this section and the next, we develop the representation theory needed for 
Grassmannian LP bounds. 

As is standard for compact Lie groups, we work with functions on Qm,n to 
find irreducible representations. Define an inner product for functions on Gm,n 
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as follows: 




f{a)g{a) da. 



Here da is the unique measure invariant on Qrn^nj normalized, so tliat J da — 1. 
That such a measure exists an is unique (the Haar measure) follows from the 
fact that Qrn,n is a compact Lie group. Equivalently, we may write 

(/,<?) := / f{U*PaU)g{U*PaU) dU, 

JU{n) 

where dU is the Haar measure on U{n), and Pa is the projection matrix for 
some fixed a G ^m,n- Now let L?'{Qm,n) denote the space of square-integrable 
functions on Qm,n- Then U{n) acts on / G Gm,n as follows: 

{Uf){Pa) ■.= f{U*PaU). 

It follows that L'^{g 

m,n) provides a representation of U(n). As we will see, this 
representation can be decomposed into irreducible subrepresentations explicity, 
and the decomposition is multiplicity-free: no irreducible representation of U{n) 
occurs more than once in L'^{Qm,n)- 

Since U{n) is a compact Lie group, its irreducible representations are well- 
studied: see for example j22l [TSj EJ [l2]. Every irreducible representation is 
indexed by a dominant weight \22\ Theorem 7.34]. In the case of U{n), we may 
take these weights to have the form [Bl Theorem 38.3] 

A = (Ai, . . . , A„) : Ai > Aa > . . . > A„, G Z. 

The dimension of the irreducible representation Vx indexed by A is given by 
Weyl's character formula [22, Theorem 7.32]. In the case of U{n), the formula 
reduces to: 

dimy.= n (3.1) 

l<i<j<n •' 

For example, the standard representation of C/(n) is indexed by A = (1, 0, . . . , 0), 
which gives 

dimV(i,o,...,o) = 

Note that there is more than one irreducible representation with the same 
dimension. 

Each dominant weight may also be thought of as a form acting on a maximal 
Abelian subgroup of the Lie group. Here A acts on the diagonal matrix d = 
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diag((ii, . . . , dn) G U (n) as follows: 



n 



1=1 



The next section describes exactly which of these forms contribute to the de- 
composition of L'^{Qm,n)- 



The group U{n)/U{m) x U{n — m) is an example of a symmetric space: a 
factor group G/K such that G is a connected semisimple Lie group and K is 
the fixed point set of an involutive automorphism of G. In this section, we use 
results from Goodman and Wallach [15] to explain how the decomposition of 
representations of Qm,n follows from this structure. 

Let Sm denote the m x m, matrix with backwards diagonal entries of 1 and 
elsewhere: 



Then U{n, s„) denotes the group of matrices which preserve the Hermitian 
form {x,y) i— > x*Sny- that is, U{n,Sn) is the set of matrices M such that 
M*SnM = Sn- This group is isomorphic the standard unitary group U{n). 
Define 



and consider the involution 9{M) := Jm,nM Jm,n on GL„(C). The fixed points 
of 9 have the form 



so the fixed point set in GL„(C) is isomorphic to GLm{C) x GLn^rn{C). 

Lemma 4.1. The Gxed point set K of 9 in G = U{n,Sn) is isomorphic to 
U{m) X U{n — m). Therefore Qm,n is a symmetric space. 

Proof. For a = (ai, . . . , am), let a denote the reversal of a, namely 



4 Symmetric spaces 






a : = 




. . . ,ai). 
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If a, b, and c have length m, n — 2m and m respectively, then we have 
Jm,n{a,b,c)'^ = {c,b,a)'^. Therefore the 1 and —1 eigenspaces of Jm,n are 
V+ = {(a, 6, a)} and V- = {(a, 0,— a)} respectively. These spaces are orthogo- 
nal with respect to the form {x,y) i-^ x*Sny- 

Now K is the set of points in U{n,Sn) which commute with Jm,n- So 
decomposing into (BV-, we have that K is the set of points in U s^) 
which leave both and V- invariant. In other words, K is the set of points 
which preserve the form Sn on the subspaces V+ and V-. Thus 

K ^U{V+,Sn\v+) X U{V-,Sn\v-) = U{n-m) x U{m). □ 

The fact that K is the fixed point set of ^ in G implies ( |151 Theorem 12.3.5]) 
that (G, K) is a spherical pair: for every irreducible representation V\ of G, the 
subspace of points fixed by K satisfies dimV^^ < 1. Those representations 
such that has dimension exactly 1 are called spherical representations. The 
following theorem [16\ Theorem V.4.3] explains how those representation relate 
to L'^{G/K). 

Theorem 4.2. Let G be a compact simply connected semisimple Lie group, 
and let K < G be the fixed point group of an involutive automorphism of G. 
Further let Gk denote the set of equivalence classes of spherical representations 
V\ of G with respect to K. Then L'^{G/K) is a multiplicity-free representation 
of G, and 

L\G/K) - V^. 

To describe which representations are spherical, we now consider diagonal 
subgroups of G and K. For d = (di, . . . ,dn), let diag(d) denote the diagonal 
matrix with diagonal entries di, . . . ,dn- Firstly, note that diag((i) is in U{n, Sn) 
if and only if dn+i-k = ^/dk-, where dk is the complex conjugate of dk- In other 
words, if d^"^ denotes the vector . . . , l/dk), then diag(d) is in C/(n, Sn) if 

and only if d = d~^. Secondly, note that if d = diag(a, 6, c) with a and c of 
length m, then 0{d) = (c, 5, a). It follows that the diagonal group 

T := {diag(ai, . . . ,am,Cm+i,- ■ .,Cn-m,am, • • • ,ai) : \ai\ = l,c = c"^} 

is contained in K. In fact, it is a maximal Abelian subgroup of K: this is called 
a torus of K. 

Recall that the irreducible representations of G are indexed by the dominant 
weights A = (Ai,...,An), where Aj > Aj+i and Aj G Z. Now the spherical 
representations of G with respect to K are indexed by those particular dominant 
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weights such that = 1 for all t = (ti, . . . ,tn) in the torus T (see Goodman 
and Wallach |15l p. 540]). So a dominant weight A is spherical if it has the 
form 

^ = (-^1) • • • 1 0, . . . , 0, —Am, . . . , — Ai) 
with Ai > . . . > Am > and Aj G Z. In other words: 

Theorem 4.3. The irreducible representations ofU{n) occurring in L^{Gm,n) 
are in one-to-one correspondence with the integer partitions with at most m 
parts. 



For any partition ^, we let H^{n), or simply H^, denote the irreducible 
representation in Lp'{Qm,n) isomorphic to o,...,o,-/i)- The Weyl character 
formula (equation (j3.ip ) now tells us the dimension of each H^. The first few 
dimensions are: 



dimF(o) = dimV(o,...,o) 
dimif(i) = dim V(i^o,... ,0,-1) 



dim i?(2) 



dinii^(i,i) 



dimF(2,i) 



dimif(fc) 



1 






1 


n^(n 


- l)(n + 3) 


n^(n 


4 

+ l)(n-3) 


(n2- 


4 

-l)2(n2-9) 
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n + k -2\'^n + 2k -I 



dimi?(i,...,i) 



n — 1 



n+l\^n-2A; + l 



/c / n + 1 



If m = 1, then ^m,n is the complex projective space CP"~^, and only the 
spaces -ff(fc) occur. In that case -ff(fe) is isomorphic to the space Harm(A;, k) of 
harmonic polynomials of homogeneous degree k in both z and z, where z = 
(zi, . . . , Zn) is a point on the unit sphere in C". Those harmonic polynomials 
were used by Delsarte, Goethals, and Seidel in their LP bounds for codes and 
designs the complex unit sphere [11 j. 

We now record a few more representations of U{n) we will need later. Given 
an nonincreasing sequence of nonnegative integers fi = (/ii,/X2, . . .), we say fi 
has size k and write |//| = A; if /i is a partition of k; that is, Yli f^i — ^- We also 
say /u has length I and write len(/i) = / if /i has I nonzero entries. For example, 



9 



(2, 1,0,.. .) has size 3 and length 2. Then for fixed Gm,n, define Hk = Hk{m, n) 
as follows: 

Hk{m,n):= H^{n). 

|/i|<fc 

len(/i)<m. 

For A; > this representation is reducible, and Hk-i is contained in Hk- When 
m = 1, Hk is isomorphic to the space of homogeneous polynomials degree k in 
both z and z on the unit sphere in C". In the next section, we will see that 
Hk is also the span of the degree-A; symmetric polynomials on the principal 
angles between a G Gm,n and some fixed h S Qm,n- Moreover, if g and h are 
polynomials in Hj^ and Hy respectively, then gh is in H/^^y, and in fact H^j^y 
is spanned by polynomials of that form. 

We also let Homfc(n) C L'^{Qm,n) denote the space of polynomials which are 
homogeneous of degree k in the entries of Pa, where Pa is the projection matrix 
of o G Gm,n- Since the constant function Pa tr(Pa) = m is in Homi(n), it 
follows that Homfc_i(n) can be embedded into Homfc(n). Similarly for fixed b, 
the distance function Pa ti{PaPb) is in Homi(n). The next section will also 
show that Hf: is a subspace of Hom^. 

James and Constantine [TT] further investigated the irreducible subspaces 
of Lp'{Qm,n)i finding zonal polynomials for each irreducible representation. We 
describe those results in the next section. 

5 Zonal polynomials 

A zonal polynomial at a point a S Qm,n is a function on points b G Qm,n which 
depends only on the the principle angles between a and h. Given any univariate 
polynomial f{x) of degree k, we define the zonal polynomial of f at b as follows: 

k 

fa{b)=Y.htr{PaP,r. 

i=0 

Here Pa and Ph are the projection matrices for a and b. As written, the zonal 
polynomial is not homogeneous, but by embedding the constant 1 into Homi (n) 
in the form tr(Pi,)/m, the exponents in fa(b) may be "pushed up" and we may 
assume fa is in Homfc(n). To see that faib) only depends on the principal angles 
between a and b, note that tr{PaPb) is simply the sum of the principal angles. 

There is another set of zonal polynomials that play a particular role in the 
theory of Delsarte bounds. Let iJ^ be an irreducible representation in Lp'{Qm,n)- 
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Then for each a S Gm,,m define the zonal orthogonal polynomial Z^^a to be the 
unique element of Hf^ such that for every p G Hf^, 

{Z^t,a,P) =p{a). 

Then zonal polynomials are invariant under the unitary group, in the following 
sense: 

Z^,b{a) = {U* Zfj^^a,U* Z^^b) = {Z^^ua, Z^^ub) = Z^^ub{Ua). 

The value of Z^^h{a) depends on the L'^(n)-orbit of (a, b) and therefore depends 
on the principle angles of a and b. With this in mind we sometimes write 
Z^^a{b) = Z^{a, b) or Z^^a{b) = . . . , ym), where (yi, . . . , Vm) are the prin- 

cipal angles of a and b. 

Schur orthogonality |22l Theorem 3.3] for irreducible representations implies 
that Z^^a and Z^^i^ are orthogonal for ^ ^ v. So, we have 

{Zfj.,a, Zu,b) = ^^i,vZ^{a,b). 

Moreover, Z^^a{b) = Z^^bict) is in fact real and symmetric in a and b. The zonal 
polynomials satisfy some other important properties, including the following 
positivity condition: 

Lemma 5.1. For any subset S C Grn,n, 

Z^{a,b)>0. 

a,b£S 

Equality holds only when XlaeS ■^m." ~ 
Proof. We have 

Z^{a,b) = Y 

a,beS a,beS 

\aes bes 
> 0. 

Equality holds if and only if J2aeS -^a*-" ~ ^- '-' 

The second important condition the zonal polynomials satisfy is called the 
addition formula: 
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Lemma 5.2. Let ei, . . . , cat be an orthonormal basis for the irreducible sub- 
space Hfj^. Then 

N 

^ei(a)ei(6) = ^^(0, 6). 
i=l 

Proof. Since Z^^a is in H^, we may write it as a linear combination of 
ei, ...,eN- 

N 
i=l 



^ei(a)ei. 



So, it follows that Z^^a(b) = ei{a)ei{b). □ 

James and Constantine give an explicit formula for the zonal orthogonal 
polynomials of Qm,n in terms of Schm' polynomials, the irreducible characters 
of SL{m, C). If y = (yi, . . . , ym) are variables and a = (si, . . . , Sm) is a partition 
into at most m parts, then the (unnormalized) Schur polynomial is defined as 



det(y--+"-^) 
det(7/. 



Each Schur polynomial is a symmetric polynomial in (yi, . . . ,ym)- For more 
information about Schur polynomials, see Stanley [23^ Chapter 7]. The nor- 
malized Schur polynomial X* is the multiple of Xq- such that X*{1, . . . , 1) = 1. 

To define the zonal orthogonal polynomials for Qm,m first define the ascend- 
ing product 

{a)s := a{a + 1) . . . (a + s - 1), 

and given a partition a = (si, . . . , Sm), define complex hypergeometric coeffi- 
cients 

m 
i=l 

Further assume we have a partial order < on partitions defined such that 
(si, . . . , Sm) ^ (ki, ■ ■ ■ , ki) if and only if Si < ki for all i. Letting y + 1 := 
(yi + 1, . . . , ym + 1), the complex hypergeometric binomial coefficients [^] are 
given by the formula 
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We can now define the zonal orthogonal polynomials for Gm,n- The following 
result is due to James and Constantine |17j . 

Theorem 5.3. Let 



Pa ■= ^Si{Si -2i + l) 



i=l 



and let a and k partition s and k respectively. Also let 



K 










a 



(k-s) 



c + 



pK—pa 

k—s 



where the summation is over partitions ai = (si, . . . , Si + 1, Sj+i, . . .) that 
are nonincreasing. Then up to normalization, the zonal orthogonal polynomial 
for is 



-x:{y), 



a<K 



where y = {yi, ■ ■ ■ ,ym) is the set of principal angles. 
The first few normalized Schur polynomials are: 
X*o{y) = 1 

^ m 

Xl{y) = -Y.yi 
m ^ 

1=1 

^1*1 (y) = -k^Yy^yj 



^ m 

XI (y) = ( + y^yj) ■ 

\ 2 ) 1=1 i<j 



Up to normalization by a constant, the first few zonal orthogonal polynomials 
are: 

Zo{y) = 1 

Zi{y) = nXl{y) - m 
Zi^i{y) = m{m - 1) - 2(n - l)(m - l)Xl{y) + (n - l)(n - 2)Xl^{y) 
Z2{y) = m{m + 1) - 2(n + l)(m + l)X*(y) + (n + l)(n + 2)X*{y). 
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A 






Absolute 
bound 




Q > 1) 


Relative 

Relative 

bound 

conditions 


n{m — a) 


n(m — a)(m — /3) 


2 

TTl — noi 

2 

m 

a < — 
n 


2 (m+l)2 (m— 1)2 / 1 1 

2{n+l) 2(n— 1) ^ ' ^' / ' rw' 

2{rn?n — Am + n) 
a + P<^ ^- ^, 

na/? m^n — 2m + n 
m"^ — 1 



Table 1: Upper bounds on \S\, when S C is an ^-code. 



The correct normahzations satisfy 

{Z^^a, Z^^a) = Zf,{l, 1, . . . , 1) = dimHf,. 

With the exception of the case /i = (which is normahzed correctly in the 
formula above), normalizations for will not play a role in the results which 
follow. 

6 Bounds 

Recall that an ^-code is a collection 5 of subspaces in Qm,n such that tr(Pa-fb) G 
A for every a ^ b in S. In this section, we find upper bounds on the size of an 
^-code in terms of either the cardinality of A or its elements. A summary of 
the results for |.A| < 2 is given in Table [TJ 

If .A = {ai, . . . , ak}, then the annihilator of A is the function 

k 

ann^(x) := - a,), 

2=1 
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The significance of tlie anniliilator is tliat ann^{tr (Pa Pb)) = for any a ^ b 
in S. More generally, for any polynomial /, an f-code is a collection S of 
subspaces such that /(tr(Pa)) ^ and f{tv{PaPb)) = for every a / 6 in S". If 
A is any set of angles and / is the annihilator of A, then an ^-code is also an 
/-code. 

Theorem 6.1. If 5 C 

Qin,n is an A-code, with \A\ — k, then 

\S\ < dim(Homfc(n)) < 

Proof. We prove more generally that if S is an /-code, with deg(/) = k, then 
l^l < dim(Homfc(n)). The result then follows by taking / to be the annihilator 
of A 

Consider the zonal polynomials fa{b) := f{tic{PaPb)), for a £ S. Note that 
fa is in Homfc(n), since /a(&) is a degree-A: polynomial in the entries of Pb- Since 
fa{b) = for every b £ S except a, and /a(a) / 0, the set {fa : a G S*} is linearly 
independent. Thus the number of functions IS"! is at most the dimension of the 
space Homfc(n). □ 



+ k-l 
k 



Corollary 6.2. Let S be a collection of subspaces in Qm,n such that tj:{PaPb) = 
a for all a ^ b in S. Then 

\S\ < n\ 

Proof. Use Theorem 16.11 with the degree- 1 annihilator of a, which induces 
zonal polynomials in Homi(n). □ 

Since fa(b) is also a degree- A; symmetric polynomial in the principal angles 
of a and b, it follows that fa is also in Hk{m,n). Then by the same argument 
as in Theorem 16. 11 we have 

Corollary 6.3. If 5 C 

Qrn,n Is an A-code, with \A\ — k, then 
\S\ < dim(IIfc(m, n)). 

If equality holds, then the functions fa form a basis for the space. Moreover, 
the space Hk{m,n) is exactly the space of functions on S. 

Theorem 16 . 1 1 and Corollarv 16.31 are called absolute bounds for Grassmannian 
codes, because the bounds depend only on the number of different inner product 
values that occur in S. When m = 1 these bounds reduce to the absolute 
bounds of Delsarte, Goethals and Seidel |1H Theorem 6.1]. There is also a 
relative bound, which depends on the actual values of the inner products and 
is sometimes tighter. 
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Theorem 6.4. Let f{yi, ■ ■ ■ , Um) G ^[ui, ■ ■ ■ , Um] be a symmetric polynomial 
such that f = '^^Cfj,Z^, where is a zonal orthogonal polynomial, and each 
Cfj, > 0. Further assume that cq is strictly positive. If S is a set of subspaces in 
Qm,n such that fa{b) := f{yi{a, 6), . . . , ym{a, b)) is nonpositive for every a ^ b 
in S, then 

Co 

Proof. Since /a(&) < for b ^ a, summing over all b £ S, we have 

Y,Ub)<fa{a) = f{l,...,l). 

bes 

Then averaging over all a G S, 

/(i,...,i)>^ Ub) 

' ' a,b&S 

By Lemma 15. H the inner sum is non- negative for fj, ^ 0. If = 0, then 
ZQ{a, b) = 1 for all a and 6, and hence, 

/(I) > 1 

= co\S\. □ 

Equality holds if and only if faib) = for every a ^ b £ S and for each 
/i / 0, we have either = or J2aeS "^M.a ~ ^- ^^^^ Section [9] that 

when > for all < deg(/), this implies that we have a Grassmannian 
t-design.) 

By way of example, we consider the case of an {a}-code in detail. 

Corollary 6.5. Let S be a subset of Qm,n such that ir{PaPb) = a for all a ^b 

in S, and a < m'^/n. Then 

n[m — a) 



\S\< 



m? — na 



Proof. The first two zonal orthogonal polynomials are ZQ{y) = 1 and (up to 
normalization) Zi(y) = X^I^Li yi—m'^/n. The annihilator for a is the polynomial 
f{x) = x — a, which induces the zonal polynomial at subspace a given by 

Ub) = tr(P,n) - = f; - a, 
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for principle angles yi, ■ ■ ■ ,ym of a and b. Thus we may write 

m 

f{yi,...,ym) = ^yi-a 



1=1 

— - a j Zo(y). 

Applying Theorem 16.41 we find that 

Co m'^/n — a 

When m = 1, we recover Delsarte, Goethals and Seidel's bound for complex 
equiangular lines: 

1 — na 

Similarly, using the zonal orthogonal polynomials Zq, Zi, Zi^i and Z2, we 
get a bound on the size of a subset containing two inner products, say a and 
p. 

Corollary 6.6. Let S be a subset of Qm,,n such that ti(PaPb) £ a, (3 for aU 
a b in S. Further assume that 

2(m^n — 4m + n) 
a + p < 



naB m?n — 2m + n 
a + l3 f < 



Then 



\S\< 



m^ Ti? — 1 

n{m — a){m — (5) 



m^ 



^ 2(n~\\ ~ + P) + 



2(ra+l) 2(n-l) 

When m = 1 this reduces to the Delsarte, Goethals and Seidel bound of 
|g| ^ n(n + l)(l-a)(l-/?) 

for lines in complex projective space CP"~^. 
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7 Other bounds 



Certain cases of equality in Corollaries 16.51 and 16.61 also achieve equality for 
bounds on the size of the largest angle in a set of subspaces. For real Grass- 
mannians, Conway, Hardin and Sloane [8] call these bounds the simplex and 
orthoplex bounds. Here we give their complex analogues. 

Recall that if Pa be the n x n projection matrix for a £ Gm,n, then Pa is 
Hermitian with trace m, so P^ = Pa — ml /n lies in a real space of dimension 
Tn? — 1. Moreover := tr(P^P^) = m(l — m/n)^ so P'^ is embedded onto 

a sphere of radius y^m(l — m/n) in M" Further recall that the chordal 
distance on Qm,n is defined by 



With this distance, the Grassmannians are isometrically embedded into ~ . 
The "Rankin bounds" given in Theorem 17.11 below (see [4^ Theorems 6.1.1 & 
6.1.2]) are bounds on the minimum distance between points on a real sphere 
as a function of the number of points and the dimension of the space. An 
equatorial simplex refers to a set of points on the unit sphere that form a 
simplex in a hyper plane of dimension — 1. 

Theorem 7.1. Given N points on a sphere of radius r in M^ , the minimum 
distance d between any two points satisfies 



Equality requires N < D + 1 and occurs if and only if the points form a regular 
equatorial simplex. For N > D + 1, the minimum distance satisfies 



and equality requires N < 2D. When N = 2D, equality occurs if and only if 
the points are the vertices of a regular orthoplex. 

Conway, Hardin and Sloane [8j apply these bounds to get the simplex and 
orthoplex bounds for real Grassmannians: we can do the same for the complex 
Grassmannians . 



dc{a,hf 



m - tv{PaPb) 





d < r\/2. 
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Corollary 7.2. Given a set S points in Gm,n, the largest inner product value 
a = {a, b) between any two points satisfies 

m\S\ — n ,„ , 

Equality requires \S\ < -n? and occurs if and only if tiie points form a regular 
equatorial simplex in R" For \S\ > r? , the largest inner product {3 satisfies 

P>—, (7.2) 
n 

and equality requires \S\ < 2(n^ — 1)- Equality occurs if the points are the 
2(n^ — 1) vertices of a regular orthoplex in 

If S is an {a}-code, then solving inequality (j7.ip for IS*! recovers the relative 
bound in Corollary 16.51 Moreover, if 1 5*1 = (equality in the absolute bound 
of Corollary 16. 2p , then 

m{mn — 1) 



a 



n 



2-1 



On the other hand, if S" is a {0, m^/nj-code, and m = n/2, then the relative 
bound in Corollary 16.61 implies that 

\S\ < 2(n2 - 1), 

which corresponds to equality in the orthoplex bound (|7.2p . 



8 Examples 

In this section we give examples demonstrating the tightness of the bounds in 
the previous sections. 

When the rank m of the Grassmannian subspaces is 1, we recover all the 
classical results of Delsarte, Goethals and Seidel ^Uj for lines in complex pro- 
jective space: their paper gives several examples of bounds with equality. In 
particular, the upper bound for {a}-codes in CP"~^ is n^, and equality can 
only hold with a trace inner product value of a = l/(n + 1). Examples of 
tightness have been found for several small values of n and are conjectured 
to exist for every n. These equiangular lines are sometimes called symmetric 
informationally complete POVMs in the quantum information literature: see 
|19j for more details or p!8] for recent results. Another important example in 
Qi^n is the relative bound (Corollary 16. 6p with inner product values of a = 
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and (3 = 1/n. The upper bound for the size of an {0, l/n}-code is n(n + 1), 
and when equahty is achieved we have what is known as a maximal set of mu- 
tually unbiased bases. Constructions achieving the bound are known when n 
is a prime power; see [H] for some constructions and |20j for apphcations to 
quantum information. 

In the case m = n/2, if (x is in Qm,n-i then its orth.ogon0,l complement d'^ is 
also in Qm^n-! and a and have a trace inner product of 0. Here again, such 
subspaces have apphcations in quantum state tomography; more details will 
be found in [13]. If is a {0,n/4}-code in Gn/2,m then by the relative bound 
(Corollarv 16. 6p . S has size at most 2(n^ — 1). In these case we may assume 
that both a and o"*" are in S, because if a and b have a trace inner product of 
n/4, then so do a-*- and b. The following construction, due to Martin Rotteler, 
demonstrates that Corollary 16.61 is tight when n is a power of 2. 

Theorem 8.1. Let Xi, . . . ,X^2_i be the Pauli matrices of order n = 2^^, and 
let 

Then U"^^^{Mj,/ — Mi} is the set of projection matrices for a {0, n/4} -code 
of size 2(n^ — 1) in Gn/2,n- 

More generally, the bound is tight when n is the order of a Hadamard 
matrix: details of the following construction will appear in in [13] . 

Theorem 8.2. Suppose there is a Hadamard matrix of order n. Then there 
exists a {0, n/4}-code of size 2(n^ — 1) in Gn/2,n- 

When the dimension of the complex space is an odd prime power, there 
is another construction which acheives the relative bound with equality. The 
following is the complex version of a set of real Grassmannian packings due to 
Calderbank, Hardin, Rains, Shor, and Sloane [7]. For lack of another reference 
in the complex case, the details are included here. 

Let V := , where q = and p is an odd prime, and let {e^ : v S F} be 
the standard basis for C". Then define the x Pauli matrices 

X{a) : Cv ^ e^+a, 
Y{a) : e„ ^ oj''^"'^ "h^ , 

where a; is a p-th primitive root of unity. Note that is an eigenvalue for 
Y{a) and e* := X^a'^*'^^" ^'^^o- eigenvalue for X{a). Define the extraspecial 
Pauli group E to be generated by all X{a),Y{a), and ljI; it has pq^ elements. 
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all of the form oj'' X {a)Y (b) , for i G Z^, a,b eV. Its center is Z{E) = {ool), and 
E := E/Z{E) is Abelian and therefore a vector space isomorphic to V'^ under 
the mapping 

(a, 6) ^ X{a)Y{h)/Z{E). 

The space V"^ has a nondegenerate alternating bilinear form (a symplectic 
form), namely 

((ai, 61), (02, &2)) := tr(af 62 - ajbi). 

It is not difficult to check that two elements in E, say w^X{ai)Y{bi) and 
X (a2)Y (b2) , commute if and only if their images in E/Z{E) satisfy 

((ai,6i), (02,62)) = 0. 

Subspaces on which the symplectic form vanishes are called totally isotropic. 
Therefore, a subspace W of E/Z{E) is totally isotropic if and only if its preim- 
age in is an Abelian subgroup. 

We now use characters of subgroups of E to define elements of Qgk^gn- Let 
W be a totally isotropic subspace of E/Z{E) of dimension n — k, and let W be 
the preimage of Vl^ in If % : — > C is a character of W , then x' : — > C 
defined by 

^{u'X{a)Y{b)) = u~\{X{a)Y{b)/Z{E)) 
is a character of W. Define a matrix 

' ' g&W 

Lemma 8.3. If W is an {n — k)-dimensional totally isotropic subspace of 

E/Z{E) and x is a character of W, then is the projection matrix for a 
q'^ -dimensional subspace of C"^" which is invariant under the action ofW. 

Proof. It is not difficult to check that 11-^ is Hcrmitian and 11^ = 11^. It is also 
not difficult to check that n^v is an eigenvector oi g ^ W for any v G C**", so 
n,^ is a projection matrix for an invariant subspace. The rank of 11,^ is the trace 
of n^, which can be computed as follows, after noting that the only elements 
of E with non-zero trace are the multiples of the identity: 

= W\i: X{9)tr{9) = -Li:--^tr(c.V) 

' ' g=ojiI i=l 

In the construction that follows we require the g-binomial coefficients, de- 
fined as 

(g" - 1) . . . - 1) 

(9--l)...(g-l) • 



n 
m 
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Theorem 8.4. For < k < n — 1, let S be the set of all q'^ -dimensional invari- 
ant subspaces of the preimages W of all {n — k) -dimensional totally isotropic 
subspaces W of E/Z{E) (as described in Lemma \K3\) . Then S is a {n — k + 1)- 



distance set in Gqk^qn of size 



,.n—k 



n 

n — k 



n (^'+1)- 

=k+l 



Proof. For j G {1,2}, let Wj be an isotropic subspace of E/Z{E), let Wj be 
its Abelian preimage in E, let Xj a character of Wj, and let Hj := 11^^. as 
in Lemma 18.31 Then 



tr(nin2 



1 



\Wi\\W2\ 

1 



Xi (51 )X2 (52)^(51 52) 



n W2I 



\Wi\\W2\ 



\Wi\\W2\ 



(or 0, depending on x'l and x'2) 



q^\WinW2\ , 
= — ^ — or . 

\Wi\\W2\ 

Furthermore, any two distinct invariant subspaces from the same isotropic Wj 
are orthogonal. If Wi / 1^2, then dim(Wi n W2) G {0, 1, . . . , n — A; — 1} and so 
I Wi n W2 1 takes n — k possible values. It follows that 5* is a (n — A; + l)-distance 
set. To find the size of S, first note that the number of isotropic subspaces of 
dimension n — fc is (see [5i Lemma 9.4.1]) 



n 

n — k 



n + 1) 

.=k+l 



and then note that each isotropic subspace produces ^ invariant subspaces. 



n 



In the case k = n — 1, Theorem 18.41 produces a 2-distance set 

of size q_i ■ The inner product values that occur are a = and /3 = : 
this construction acheives equality in the relative bound (Corollary 16. 6p . In 
his thesis, Zauner ^25j has a construction which has these same parameters (in 
fact, Zauner 's construction is more general, as it also allows q to be an even 
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prime power). In the case fc = n — 2, we get a 3-distance set in Qg„~2 gn of size 

(g^-i)(g-i) — ^ with inner product values a = 0, P = q"" , and j = q'^ ■ 
There are many open questions regarding whether or not tightness in the 

bounds can be achieved; in particular, it is not known if there are any examples 

of subspaces achieving equality in the absolute bound (Corollary [62]) for m > 1. 

The smallest nontrivial case is a set of 16 subspaces of dimension 2 in C^, with 

an inner product value of a = 14/15. 



9 Designs 

In this section, we introduce the concept of a complex Grassmannian 2-design. 
We give lower bounds for the size of a t-design and indicate the relationship 
between designs and codes. 

Recall that Ht{m, n) is the direct sum of the irreducible representations 
of U{n) containing the zonal polynomials Z^^a, where fi is an integer partition 
of size at most t and length at most m. Ht{m, n) may also be thought of as the 
symmetric polynomials of degree at most t in the principle angles of pairs of 
subspaces in Qrn^ri' Since the zonal orthogonal polynomials Z^^^ (with ^ t 
and len(^) < m) span Ht{m,n) and are contained in Homt(n), it follows that 
Ht{m,n) is a subspace of Homt(n). 

We call a finite subset S C Qm,n a t-design if, for every polynomial / in 
Ht{m,n), 




/(c) dc. 



In other words, the average of / over S is the same as the average of / over 
the entire Grassmannian space. Recall that the average of / over Gm,n can be 
written as (1, /): with this in mind we define an inner product for functions on 
S as follows: ^ 

if'9)s ■= ■r^y"/(«)5(o)- 

Then S is a t-design if (1,/) = {l,f)g for every / G Ht{m,n). Equivalently, 
the zonal orthogonal polynomials Z^^a span H^, so S* is a t-design if every Z^^a 
has the same averages over S and Gm,n, where /i is a partition of at most t into 
at most m parts. 

By way of example, consider Theorem 16.41 If / = c^g^ and > for 
every \fi\ <t, then equality in Theorem 16.41 implies that S" is a t-design. 
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For the purposes of quantum tomography apphcations, 1- and 2-designs 
play a special role (see [13], as well as [21]). In those cases, there is a more 
explicit description of a i-design. 

Lemma 9.1. Let S be a finite subset of Gm,n- Then S is a 1-design if and only 
if 



Moreover, S is a 2-design if and only if 



1-^2 Pa® Pa = I Pa® Pa da. (9.1) 



1^1 ... 

Before proving Lemma [9. H we note the integral on the RHS of equation (j9.ip 
can be evaluated explicitly. Writing Pa = YliLi o-iO-l for some orthonormal basis 
{oj} of a, and letting T denote the "swap" operator T : e-i ® Cj ^ ej ej, the 
integral is obtained from Lemma 5.3 of |20j : 



Pa ®Pada = — ^ [{nm - 1)1 + (n - m)T] . 

nyw^ — 1) 

Proof. We prove the lemma by showing that Ht = Hom^ for t G {1,2}; the 
result then follows by considering the polynomials of the form a i— > {Pa)ij in 
Homi and a i-^ {Pa)ij{Pa)ki in Hom2. 

Recall that Ht is contained in Homf, so it suffices to show that the dimen- 
sions of the spaces are equal. When t = 1, we have dim(i/i) = dim(Homi) = n^, 
so Hi = Homi. When t = 2, recall that dim(//2) = C2 ) (assuming m > 1), 
and the space of homogeneous degree-2 polynomials on the coordinates of n x n 
matrices has dimension . However, Hom2 is the space of degree-2 polyno- 

mials on projection matrices, not general matrices. If Pa is a projection matrix, 
then the degree-2 polynomial 

^ mtr{AP^) - tr(P,) tr(AP,) 

is identically zero for every A. There are linearly independent polynomials 
of that form for general n x n matrices; therefore, 

J- /TT N fn^ + '^\ 2 
dim(Hom2j =1 1 — n = 



2 

Thus H2 = Hom2. □ 
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We now consider bounds for t-designs. The following is the absolute bound. 
Lemma 9.2. If S is a t-design, then 

\S\ > dim(i?Lt/2j(m,n)). 

Proof. Let {ei,...,eAr} be an orthonormal basis for Since Cj is a 

symmetric polynomial in the eigenvalues, so is ejcj. It follows from the unique 
decomposition of L?'{Gm,n) that ejcj is in H2\t/2\ ^■iid therefore in H^. If 5 is a 
t-design, and ejcj is in Ht, then 

{ei,ej) = (IjCjej) = {l,eiej)g = {ei,ej)g, 

whence it follows that {ei, . . . , e;} are orthogonal as functions of S (a space of 
dimension 151). □ 

If equality holds, then the basis for Ht/2{'m, n) is also a basis for the functions 
on S. There is also a relative bound. 

Theorem 9.3. Let f{xi, . . . ,Xm) G '^[x] be a symmetric polynomial such 
that f = ^^c^Zfj,, where is a zonal polynomial for the Grassmanian 

space, and cq > 0. Furthermore, suppose S is a t-design such that fa{b) = 
f{yi{a, b), . . . , ym{0', b)) > for every a ^ b in S, and C/^ < for every > t. 
Then 

|S| > ffiiiiiii). 

Co 

Proof. Let fa be the zonal polynomial of / at a, so that fa{b) > for 6 7^ a. 
Summing over all b & S, 

\S\{lJa)s>fa{a) = f{l,...,l). 

Again averaging over all a G 5, 

fil,...,l)<^{l,fa)s 

a&S 

aeS M 
/u aeS 
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Since S is a t-design, the inner sum is zero for |/i| <t (fi ^ 0). For |^| > t, the 
inner sum is nonnegative (by Lemma IS.ip and < 0. Therefore, 

aeS 

= co|5|. □ 

If equaUty holds, then we have /(a, 6) = for every a ^ b in S. That 
is, S is an /-code. Furthermore, for every > t, we have either = or 

As with classical codes and designs, the case where S is both a /-code and 
a t-design is of particular interest, as the size of the set can be determined 
exactly. Combining Theorems 16.41 and 19.31 gives the following. 

Theorem 9.4. Suppose S is an f-code for f = c^Z^, where > 0, and S 
is also a t-design for t > deg(/). Then 

I5|_/(M.....1) 



Co 

Consider the following polynomial in Ht{m,n): 



len{fi)<m 

This polynomial satisfies {Zt^a, f) = f{o) for every / G Ht{m, n). Taking f = Zt 
in Theorem 19.41 we get: 

Corollary 9.5. If S is a Zf-code and a 2t-design, then 

\S\ = dini(H-t{m,n)). 

Theorem 9.6. Any two of the foUowing imply the third: 

• S is an f-code, where deg(/) = t; 

• S is a 2t-design; 

• \S\ = dim{Ht{m,n)). 

Proof. Suppose 5 is a /-code with \S\ = dim{Ht). Since equality holds in 
Corollary [631 the polynomials fa are a basis for Ht. However, we have 



{Zt,a,fb) = fh{a) 



lO, 6 /a; 

\/(l,l,...,l), b = a. 
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Thus {Zf^a} is a dual basis for Ht and each Zt^a is a multiple of ft^a- Now 
consider the averages {Zt^a, fb)s- since /a(&) = Zt^aip) = for 6 7^ a, we get 



(■^t,a, fb)s 

Thus we have 



0, 6 / a; 

/(1,1,...,1), 6 = a. 



(1, Zt^afb)s = (^i,ai /fe/g — \-^t,a, fb) — (l, Zt^afb) 

for the bases {Zt^a} and {/b}. But the set {Zf^afb} spans H2t{n), so S" is a 
2i-design. 

Conversely, suppose 5 is a 2t-design with \S\ = dim{Ht), and let / annihilate 
of the angle set of A. Since Ht spans the functions on IS*!, each fa is in Ht and 
is therefore a polynomial of degree t. Thus / has degree t. □ 

The simplest case of Theorem 19.61 is when t = 1: in this case, is a 1- 
distance set and a 2-design of size n^. Moreover, S" is a Zi-code, and Zi is 
the annihilator of "^^^""^-^ . Thus the inner product between every two distinct 

1 • m(mn—l) 

subspaces is q = ■ 



10 Association schemes 

As Theorem I9.6l indicates. sets of Grassmannian subspaces which reach equality 
in the Delsarte bounds have a great deal of structure. In this section, we show 
that — much like spherical codes and spherical designs — these sets are often 
endowed with the structure of an association scheme. 

Let S be an /-code with a finite number of distinct sets of principal angles 
y = (yi, . . . , Urn)- Denote the set of y's that occur by y. For each y £ y, define 
a \S\ X IS*! matrix as follows: 

I 1, a, 6 have principal angles y; 
Ay{a, b) := < 

I 0, otherwise. 

Each Ay is a symmetric {0, l}-matrix. Furthermore, each pair (a, b) has some 
principal angle y, so J^yey^y ~ ^ ■• '^h^re J is the all-ones matrix. If := 
(1, . . . , 1) denotes the trivial principal angles set, then Aq := Ay^ is the identity 
matrix. We will call the Ay matrices Schur idempotents, as they are idempotent 
under Schur multiplication, defined as follows: 

{AoB),j := AijBij. 
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Under certain conditions, these Schur idempotents form an association scheme. 

For each integer partition ^ and corresponding zonal polynomial Z^, define 
an |5| X \S\ matrix as follows: 

Each En is also symmetric and in the span of {Ay}y^y: 

' ' yey 

In particular, Eo is a scalar multiple of J. When {Ay}y^y forms an association 
scheme, the matrices E^ are the scheme's idempotents. 

Lemma 10.1. If S is a 2t-design, then {-E//}|y:/|<t,ien(;i)<m ^.re a set of orthogonal 
idempotents. 

Proof. Suppose = i and |A| = j, with i,j < t. Then 

{Ei,Ex)a,b = J^Y1 ^M^"' c)^a(c, b) 

= 1^ {^tJ.,a, Zx,b)s ■ 

Since Z^^a and 5 arc in Hf, their product is in H2t- Now is a 2t-design, so 
the average of Z^^aZx^b over S is the same as the average over Qm,n- But 

{Zij.,a,Zx,b) = Sx,i^Zf,{a,b), 
and so we find that Efj,Ex = ^x^^E^. □ 

More generally, if = i and |A| = j, and 5 is a (i + j)-design, then E^ and 
Ex are orthogonal. 

Now suppose is a 2t-design. By the previous lemma {£'yii}|^|<i are linearly 
independent, and clearly the matrices {Ay}y(zy are also linearly independent. If 
|3^| equals the number of partitions of at most t (into at most m parts), then the 
span of {Ay^y^y and {£^^}|^|<i are the same. Since {£^^}|^|<4 is closed under 
multiplication, so too is the span of {Ay}y^y, and so we have an association 
scheme. 

Corollary 10.2. Let S be a 2t-design in Qra^n with principal angle set 3^. If 
\y\ is equal to the total number of partitions of 0,1, ... ,t into at most m parts, 
then {Ay}yQy is an association scheme. 
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Lemma 10.3. Let S be a 2t-design in Qm^n with principcil angle set 3^ such 
that \y\ is the total number of partitions of 0,1, ... ,t into at most m parts. 
Then {-E^}|^|<t,ien(^t)<m the idempotents of the scheme {Ay}y^y. 

Proof. Since = j^Yly^y Zt^{y)^y' see that the matrix [Z^{y)] is 
the transition matrix between the two bases of the association scheme and is 
therefore invertible. It fohows that for each yi in 3^, some Hnear combination of 
the rows forms a homogeneous degree-t polynomial Qi such that gi{yj) = 5ij. 
(Conversely, if such gi polynomials exist, then [Z^[y)] is invertible.) Then 

c:y{a,c)=yi 
= {9i,a, Zf^fi) ■ 

Now write gi = I]|A|<t Ci,\Zx, so that 

\\\<t 

Thus AiE^ = Ci^^E^ for some q,^. □ 

By way of example, let t = 1, and suppose S" is a 2-design with only one 
nontrivial principal angle set (and one trivial one, for a total of two). The 
number of partitions of at most 1 is also two (/x = and = (1)), so by 
Corollary 110.21 we have an association scheme. In this case the scheme is the 
trivial one, namely {/, J — I}. 

As another example of an association scheme obtained from principal angles, 
consider the collection of subspaces in Gn/2,n from Theorem l8.ll This collection 
has four distinct sets of principal angles: 

y = (1, . . . , 1) (trivial principal angles), 
y = (0, . . . , 0) (angles between a and a"*"), 
2/= (1^_^,0^_^), 

n/4 n/4 

y = (i...,^). 

While |3^| = 4 is the number of partitions of at most 2 (/x = 0, = (1), 
fi = (1,1) and fi = (2)), the hypotheses of Corollary 110.21 are not satisfied 
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because the subspaces do not form a 4-design. Nevertheless, it is easy to verify 
computationally that this collection does give a 3-class association scheme. 

We may define a coarser set of relations on an /-code S using the sums of 
principal angles — the inner products of the projection matrices — instead of the 
principal angles themselves. Let A denote the set of nontrivial inner product 
values that occur in S, so S is an ^-code. For a G ^ let A'^ be the \S\ x \S\ 
matrix defined as follows: 

A'{ab)-=l^' tr(P„Pfe)=a, 
" 1 0, otherwise. 

Also define A'^ := I for the identity relation. Clearly each A'^ is in the span of 
{Ay : y £ y}; in fact 

In particular, A'^ = Aq = I, and if is in A, then Aq = ^(o,...,o)- before, the 
matrices are Schur idempotents and sum to J. Next we need the corresponding 
idempotents. For each i £ {0, . . . ,t}, define E'- as follows: 

|/i|=i 

This implies that E'q = J/\S\ and E'-{a,b) = {Zi{a,b) — Zi-i{a,b))/\S\ for 
i > 0. As in Lemma 110.11 if S" is a 2t-design, then {E'- : i < t} is a set of 
orthogonal idempotents, and if is a {2t — l)-design, then {E[ : i < t} are 
linearly independent. 

Clearly E'^^ is in the span of {Ay : y G 3^}, since each E^ is in that span. But 
suppose Zi{y) is the annihilator polynomial of some i-distance set, so it is a only 
function of yi: then in fact E'^' is in the span of {A'^ : a £ A}. If Zi{y) is an 
annihilator for sufficiently many i, then {£"• : < i <t} and {A'^ : a £ Au{m}} 
span the same set, and that set is closed under multiplication. 

Corollary 10.4. Let S be a 2t-design that is also an A-code in Qm,n- 

and Zi{y) is an annihilator polynomial for each i <t, then {A'^ : a £ AU {m}} 

is an association scheme. 

In fact, these hypotheses can be weakened. 

Theorem 10.5. Let S be a {2t — 2)-design that is also an A-code in Gm,n- If 
\A\ = t, and Zi[y) is an annihilator for each < i < t — 1, then {A'^ : a £ 
A U {?Ti}} is an association scheme. 
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Proof. Since 5 is a 2{t — l)-design, the idempotents {E'^ : < i < t — 1} 
are linearly independent. We claim that / is also linearly independent from 
{E'- : < i < t — 1}. For, if / = Yll=o CiE[^ then the off-diagonal entries of / 



But all off-diagonal entries are 0, implying that the polynomial has t roots in 
Ui, a contradition. So {E[ :0<z<t — 1}U {/} is linearly independent and 
therefore spans {A'^ : a G {m}}. Since it is closed under multiplication, we 



By way of example, suppose t = 2 in Theorem I10.5[ Note that ^o(y) and 
Zi{y) are always annihilators. It follows that if 5 is a 2-design, and the inner 
product set A = {tr(i-*aP{,) : a 7^ 6 G 5} contains exactly two distinct values, 
then {A^ : a G ^ U is a 2-class association scheme. 

Corollary 10.6. Let S he a {2t — 2)-design and an A-code in Qm,n such that 
\A\ = t and Zi[y) is an annihilator for i < t — 1. Then the idempotents of the 
scheme {A'^ : a £ AU {m}} are E'q, . . . , E[_i, and J — YllZo E[. 

Proof. Let fa denote the annihilator polynomial of A\{aQ,a}, normalized 
so that fa{(y) = 1. Then is a polynomial of degree t — 1 in Vi^ a-iid the 
corresponding zonal polynomial f^^a is in Ht-i{n). Writing Pi := Zj — = 
E|m|=» we have 




have an association scheme. 



□ 




tr(PaPc)=OL 



= {fa,a,Pi,b)s - ^^^nMPaPb)) 
= {fa,a,Pi,b) - ^4^mtr{PaPb))- 





= Ca,iPi{tT{PaPb)) - ^^^J^P^{tv{PaPb)) 
= {Ca,i\S\ - fa{m)){E'j)a,b- 



Thus A'^E[ = Xa,iE[ for some constant A, 



□ 
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